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Abstract 

We compute the effects of soft gluon resummation for the production of 
high mass systems in hadronic collisions. We carefully analyse the growth of the 
perturbative expansion coefficients of the resummation formula. We propose an 
expression consistent with the known leading and next-to-leading resummation 
results, in which the coefficients grow much less than factorially. We apply 
our formula to Drell-Yan pair production, heavy flavour production, and the 
production of high invariant mass jet pairs in hadronic collisions. We find that, 
with our formula, resummation effects become important only fairly close to 
the threshold region. In the case of heavy flavour production we find that 
resummation effects are small in the experimental configurations of practical 
interest. 
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1 Introduction 



In this work we deal with the problem of the resummation of logarithmically 
enhanced effects in the vicinity of the threshold region in hard hadroproduction pro- 
cesses. Drell-Yan lepton pair production has been in the past the best studied ex- 
ample of this sort 0, [|. The threshold region is reached when the invariant mass 
of the lepton pair approaches the total available energy. A large amount of theoret- 
ical and phenomenological work has been done on this subject. References Q and 
H summarize all the theoretical progress performed in this field. Resummation for- 
mulae have also been used in estimating heavy flavour production @, 0. In this 
case only a leading logarithmic resummation formula is known. Calculations of the 
next-to- leading logarithms are in progress ||. 

In the present work, we will mostly be concerned with difficulties that arise when 
one tries to apply resummed formulae to physical processes. This is a highly non- 
trivial problem. What one typically finds is that resummation involves the integration 
of the running coupling over the Landau pole, which has to be regulated. In early 
works on resummation in the Drell-Yan cross section, the problem was avoided by 
either assuming a fixed coupling constant (Curci and Greco in ref. ]]J), or by shifting 
the argument of as so as to move the position of the Landau pole to Q 2 = (Parisi 
and Petronzio in ref. Jl|). In ref. || a cut-off procedure was introduced in order to 
regulate these singularities. A similar approach was used in ref. ||, in the context of 
heavy flavour hadroproduction. In refs. pi and |Tl| a principal value prescription 
was adopted instead, and an application to top production was proposed in ref. 0. 
It was generally found that threshold corrections become quite large, long before the 
hadronic threshold region is reached. 

Lately, the problem of the presence of an integration over the Landau pole in re- 
summation formulae has been reexamined from the point of view of the occurrence of 
infrared renormalons (IR) in the QCD perturbative expansion for the Drell-Yan pro- 
cess [0, [12|, 13]. Roughly speaking, the IR point of view relates the factorial growth 



of the coefficients of the perturbative expansion to the presence of power-suppressed 
corrections to the process in question. It is found that the ambiguities associated 
to the resummation of the non- convergent (asymptotic) perturbative expansion have 
precisely the form of a power-suppressed correction. In this context the separation 
of perturbative and non-perturbative effects is at best ambiguous, since it relies on a 



specific prescription needed to resum an asymptotic expansion. In refs. |L0|, [T^, [13 1 
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it was argued that the Landau pole integration in the resummation formulae for 
threshold corrections gives the leading IR behaviour of the perturbative expansion 
of the Drell-Yan cross section, and that the associated factorial growth is the one 



corresponding to a 1/Q ambiguity. A recent work of Beneke and Braun [13] has 
however demonstrated that the approximations made in the resummation formulae 
for the logarithmic corrections are insufficient to correctly describe the IR structure, 
and that when higher order contributions are properly included no factorially growing 
terms associated to a 1/Q ambiguity do arise in the Drell-Yan cross section. This 



result was subsequently confirmed in refs. [Ijj and fl6fl . In the latter reference, the 
absence of 1/Q effects was shown to be a consequence of cancellations related to the 
Kinoshita-Lee-Nauenberg theorem. 

The IR point of view teaches us a very important lesson to keep in mind: it 
is not enough to make sure that all the leading corrections are properly included 
in the perturbative expansion. We should also make sure that formally subleading 
terms, which are not controlled in our approximation, will not affect the asymptotic 
property of the expansion. In fact, if formally subleading terms happen to have a 
strong factorial growth, they may induce large corrections even in kinematic regions 
where our resummation is not justified. 

In the present work we find that besides the IR problem, other, more important 
sources of factorial growth may be introduced, which are spurious and are by no means 
implied by the threshold approximation, since they are not enhanced at threshold. 
These large terms arise when one attempts to formulate the resummation problem 
in x-space, as opposed to its natural formulation in moment (or N) space. To be 
specific, let us consider the case of the Drell-Yan pair production. The Drell-Yan 
cross section can be written in x space (schematically) as 

<r( DY )(r) = f 1 dxt dx 2 dx F(xi) F(x 2 ) <r (DY) (x)S(x 1 x 2 x - r) . (1.1) 
Jo 

In moment space a simple factorized expression follows 



(DY) _ „ p ~(DY) n 9 x 



where 



1 ^ X N F{x) 
x 



JO T 
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1 dx 



a^ BY \x) . 



(1.3) 



f — x 

X 



The threshold region r — > 1 corresponds to the limit iV — > oo in A^-space. In this 
limit, soft-gluon radiation produces large logarithmic corrections of the type Og ln m A" 
that are resummed in the partonic cross section a. Resummation of soft-gluon effects 
is best expressed in iV-moment space, because it leads to the exponentiation of the 
logarithmic corrections. Exponentiation is a consequence of dynamics and kinematics 
factorization. By dynamics we mean factorization of multigluon QCD amplitudes to 
logarithmic accuracy. By kinematics we mean factorization of the phase space: the 
constraint of longitudinal-momentum conservation factorizes in iV-moment space. 

However, the moment space formula can be turned to an rr-space formula. We will 
see that with this transformation, by neglecting certain subleading terms, one may 
generate large, factorially growing corrections, which may be wrongfully attributed 
to the original resummation formula. In fact these subleading terms are there to 
compensate for the fact that exponentiation is imperfect in the x-space formulation, 
and should not be dropped. If they are neglected, kinematic constraints that were 
satisfied in the original formulae (e.g. momentum sum rules) are strongly violated, 
and actually diverge, in the x-space expression. A typical consequence of this proce- 
dure is that the final formula for the physical cross section receives large soft gluon 
corrections (actually, divergent ones) even if we are far from the threshold region. 
With these large factorially growing terms, a corresponding power corrections of the 
order of (Aqqd/Q) 5 can be associated, where 5 can be much less than I. The usual 
"common sense" assumption that power corrections arise from regions of phase space 
where the momenta are of order Aqcd seems therefore to fail, and one is forced to use 
a cut-off of several GeV in order to make any sense out of the resulting formulae. We 
argue that all these paradoxes are simply avoided if the transformation to x-space is 
performed exactly. 

In the present work we provide a specific prescription for the resummation of soft 
gluon effects that does not have any factorially growing terms in its perturbative 
expansion. The ambiguity associated with the perturbative expansion of our formula 
is therefore free of any 1/Q M effects for any M > 1, and in fact turns out to be of 
the form e~ H ^^^ T ^ A , where r is the ratio of the squared invariant mass Q 2 of the 
produced heavy system and of the total CM energy squared, H is a slowly varying 
positive function, and A is the QCD scale parameter. 

The fact that the ambiguities introduced by the Landau pole only arise for values of 
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1 — t of the order of A/Q is consistent with our expectation that a correct resummation 
of soft logarithms should allow control over the perturbative expansion down to the 
scale at which the coupling constant blows up, namely A. A successful resummation 
program, in fact, should work regardless of the size of the logarithmic terms being 
resummed, provided one can prove that the neglected terms are sufficiently small. 
The non-perturbative regime is not defined by the region of momenta in which higher 
order terms are larger than lower order ones; it is defined by the domain in which 
the terms which are neglected by the resummation procedure are comparable in size 
with those taken into account. The universality of soft gluon emission should allow 
full control to be maintained over the dominant contributions to the perturbative 
expansion even when they become large, and to correctly resum them regardless of 
their size. Our approach to resummation shows that this is indeed possible, confining 
the effects of the really non-perturbative regime of QCD to their natural scale, namely 
A. 

The paper is organized as follows. In section 2 we will give a few reference for- 
mulae and establish our notation. We will mainly deal with the Drell-Yan case as 
an illustrative example. In section 3 we will show how large spurious corrections may 
arise in the computation of the resummed Drell-Yan cross sections. In section 4 we 
will propose an alternative resummation method, and in section 5 we discuss its im- 
plications for Drell-Yan cross sections. In sections 6 and 7 we discuss heavy flavour 
production, and in section 8 jet production at large transverse momenta. In section 9 
we discuss a few remaining issues, and in section 10 we give our conclusions. In 
Appendix A we prove the absence of factorially growing terms in our resummation 
formula, in Appendix B we discuss some details about the numerical method we used, 
and in Appendix C we derive some results regarding the inverse Mellin transform in 
the leading logarithmic approximation. 

2 Basic formulae and notation 

We begin with the formula for Drell-Yan pair production 

a(T, Q 2 ) = f 1 dx dx! dx 2 F(x 1 )F(x 2 )5(xx l X2 - t)A(x, Q 2 ) . (2.1) 
Jo 

At the Born level, omitting obvious factors, we have A(x,Q 2 ) = 5(1 — x). We also 
omit, for ease of notation, the parton indices. The cross section is given in terms of 
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the parton densities F(x) as measured in the deep inelastic processes. Denning the 
Mellin transform as 

f N =f 1 -z N f(z) (2.2) 
Jo z 



we can rewrite eq. (|2.1|) as 



a N (Q 2 ) = F^(Q 2 ) A N (Q 2 ) . 



(2.3) 



The resummed coefficient function for the Drell-Yan process in the DIS scheme is 

SSI 



lnA 7V (Q 2 



dx 



x 



N 



1 — X 



2 W„2 



(l-x)Q 
(l-:r) 2 Q 2 g 2 



-A(a s (q 2 )) + B(a s ((l-x)Q 2 )) 



0(a s (a s \nNY 



with 



7T 



A(a s ) = — + ( ^) A® , B (a s ) = ^flW 
where [G A = 3 , C F = 4/3 , T R = 1/2 in QCD) 

A® = C F , 
and the coefficient if is given by J17 



A® = -C F K , flW = 

2 ' 2 



67 7T 



10. 



K = C A \—- — \- —T R N f . 



18 6 



9 



(2.4) 
(2.5) 

(2.6) 
(2.7) 



Note that, due to the integration of the running coupling, the integral in eq. (|2.4j ) 
is singular for all values of N. However, if we perform the integration to next-to- 
leading logarithmic (NLL) accuracy (i.e. we compute the leading a 7 ^ ln n+1 N and 
next-to- leading ln n N terms) [|| , we find 

\nA N (Q 2 ) = IniV gi (b a s \nN) + g 2 (b a s \nN) + 0(a k s ln fe-1 N) (2.8) 

where a$ = «s(Q 2 ), &o an d &i are the first two coefficients of the QCD /3-function 



UC A -4T R N f , 17C 2 A -10C A T R N f -QC F T R N f 



(2.9) 



12tt " 1 24tt 2 

and the leading and next-to-leading functions g\ and g 2 are given by (je — 0.5772 . . . 
is the Euler constant) 



irb \ 



[l - 2A) ln(l - 2A) - 2(1 - A) ln(l - A) 



(2.10) 



-6- 



52(A) 



+ 



+ 



7T 2 6§ 
7T&jj 



21n(l - A) - ln(l - 2A) 



ln(l - A) + 



ln(l - 2A) -21n(l 



ln(l - A) - ln(l - 2A) 
2A) 



(2.11) 



A) + -W(l 



ln 2 (l 



A) 



We observe that the A^-space formula eq. ( |2.8| ) is finite and uniquely defined up to the 
very large value of N = Nl = exp 



2a$bo 



in spite of the fact that it is obtained from 
an expression which is formally divergent for any value of N. It was shown in ref. || 
that the expression x N — 1 in formula ( |2.4j ) can be replaced to NLL accuracy by the 
theta function —9(1 — x — e~ lE /N). The region where the integral is divergent is 
therefore excluded from the integration if AT < Nl. This shows that the divergences 
present in the integral of eq. (|2.4j) are subleading for large N. They may be cancelled 
by other divergences of the same nature, neglected by the approximations that lead 
to formula (|2.4j ). This is indeed the case for the leading 1/Q singularity, as shown in 



ref. [p 



Our result in eq. Q2.8|) can be expressed in terms of as(/z 2 ) for an arbitrary value 
of the renormalization scale fi 2 . We thus achieve full control over the renormalization- 
scale dependence. To do that, we must take into account the scale dependence of the 
next-to- leading function g 2 , which is given by 



where g2(\,Q 2 ) is (72(A) as defined above and 



Si (A) 



TT&qA 2 



ln(l - 2A) - 21n(l - A) 



(2.12) 



(2.13) 



Also the factorization- scheme dependence is completely under control up to NLL 
accuracy fL8|l . In the MS scheme the resummed coefficient function can still be ex- 
panded as in Eq. fl2.8|) , that is: 



lnA5f(Q 2 ) = lniV^ s (6 oa5 lniV) +^ 2 MS (6 oa5 lniV) + 0(4 h/" 1 AO , (2.14) 
and the leading and next-to-leading functions g^ 13 and g^ s are given by 



<?f S (A) 



AW 

TlbnX 



2A + (l-2A)ln(l-2A) 



(2.15) 
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^ 2 MS (A) 



A (2) 



7T 



2A + ln(l - 2A) 



2A0) 



IE 



Tlbn 



ln(l-2A) 



+ 



nb 3 



2A + ln(l - 2A) + - ln 2 (l - 2A) 



(2.16) 



In this scheme, the renormalization-scale dependence is again taken into account by 
the function g 2 as in eq. ( |2.12| ), and the corresponding function g[ is: 



<7i (A) 



vrfenA 2 



ln(l - 2A) + 2A 



(2.17) 



Sometimes, as an illustration, we will also use the double log approximations 
(DLA) to the resummation formulae 

A 0) 



<?i(A) 



Si MS (A) 

which give rise to the expressions 

AO) 

\nA N (Q 2 ) 
lnAf^Q 2 ) 



6 7T 

2 AO) 



\ + 0{\ 2 ) , 
A + 0(A 2 ) , 



7T 

2 AO) 



7T 



7T&t 

a s In 2 N + 0((a s In A^) fc+1 In N) , 
a s In 2 N + 0{{a s In A^) fc+1 In N) 



(2.18) 
(2.19) 

(2.20) 
(2.21) 



3 Problems with x-space resummation formulae 

In this section we discuss the problems that may arise when turning the iV-space 
resummation formula into an x-space formula. We begin by considering the DLA case. 
No running coupling effect is present in this case. Therefore, there is no Landau pole, 
and the result of resummation should be finite and free of ambiguities. There is a 
"realistic" limit that corresponds to this case, which is the limit of large colour factors 
and small coupling. In the case of gluon initiated processes in the MS scheme, this 
limit is not far from reality. For the sake of definiteness, we also fix the structure 
function 

6 

Fn = N(N + l)(N + 2) (3,1) 
which corresponds to the x-space structure function 

F(x) = (l-x) 2 . (3.2) 
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Our cross section formula is 

1 rC'+ioo 

<t{t) = — F 2 N A N r- N dN, (3.3) 

2m Jc-ioo 

and we now take the simplified soft gluon factor to be given by 

A N = exp(alog 2 iV) , (3.4) 

where only the double logarithmic term has been kept in Ajy. The coefficient a is 
both process and scheme dependent (see Sects. |] and §); it is given by a = CpO-sl^ 
and a = 2Cfois/^ for the qq initial state in the DIS and MS scheme respectively, and 
by a = 2CaoisI^ for the gg initial state in the MS scheme. The integral ( |3.3| ) is not 
absolutely convergent for large N, since An grows faster than any power for large N. 
Observe, however, that if we expand A at in powers of a, the integral converges order 
by order in perturbation theory. In the corresponding perturbative expansion, we can 
therefore deform the integration contour into two straight half-lines from C — (i + e) oo 
to C, and then to C + (i — e)oo. Once the perturbative expansion is written in this 
way we realize that it can be resummed into the expression 

1 r G+(i-e)oo 

a(r) = — F*A N r- N dN, (3.5) 

2m Jc-(i+e)oo 

which is now convergent, and independent of e. In the following, we will always 
interpret the integral in the above sense, omitting the explicit reference to the e. 

Let us now see what happens if we try to rewrite eq. ( |3.3| ) as an x-space formula. 
First of all we need the inverse Mellin transform of An- We have (see Appendix C 
and ref. 11) 

"j rC+ioo 

A(x) = / exp(alog 2 AO x~ N dN 

2m Jc-ioo 

= (0(1 -rj-x) exp[alog 2 (l - x)]\ x (1 + NLL terms) (3.6) 

where, according to the usual definitions, NLL terms stand for contributions of the 
form a k s log" 1 1/(1 — x) with k > 1 and m < k. The role of the 9 function is to define 
the right-hand side in a distribution sense as rj — > 0, so that we have the correct 
normalization of the first moment 

l 

A(x) dx = 1 . (3.7) 



o 
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Thus, defining 

C(z) = f^F (A) F (x 2 ) (3.8) 

Jz X 2 \X 2 J 

and neglecting NLL terms, we obtain the following expression for the cross section 

a ( T ) = / dxdx\ dx 2 F(xi) F(x 2 )S(xxi x 2 — r) A(x) 
Jo 

= C dx exp[alog 2 (l - x)]—C (-) , (3.9) 
J t dx \ x J 



where an integration by part was performed. We now see that the integral in eq. ( |3.9| ) 
is divergent at x = 1 for any value of r, since the expression exp[a log 2 (l — x)] diverges 
faster than any power as x — > 1. Let us examine more carefully the origin of this 
divergence. If we expand eq. ( |3.9| ) in powers of a, each term of the expansion is 
integrable, but the corresponding series is divergent. Since the term d/dx C{r/x) is a 
smooth function of x as x — > 1, the nature of the divergence is given by the following 
integral 

/ 1 exp[alog 2 (l-x)]dx = f; £ C log 2 ^ dz = £ ^11. (3.10) 

JO k=Q K\ Jo k=Q Kl 

The asymptotic behaviour of the coefficients for large k is (2k)\/k\ ~ 4 k k\. The 
expansion is therefore an asymptotic one. Observe that the lower limit of the integral 
in eq. (|3.10|) is irrelevant for this conclusion. It is known that factorially growing 
terms in the perturbative expansion are associated to power-like ambiguities in the 
resummed expression. In order to resum the asymptotic expansion, we should in fact 
truncate the series when the next term is of the same size as the current one, i.e. 
when 4ak = 1. The error on the resummed expression is then of the order of the left 
over term 

5= (Aa) k k\ w (Aa) k k k e- k = . (3.11) 
If we replace the appropriate value of a we get 

5 = (-) qq, DIS, (3.12) 

qq, MS, (3.13) 
99, MS. (3.14) 
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Although it is power-suppressed, the smallness of the exponent makes this effect 
potentially large. Thus, for example, for rif = 5 we have for Drell-Yan in the DIS 
scheme a (A/Q) ' 72 power correction, and for heavy flavour production via gluon fusion 
at fixed invariant mass of the heavy flavoured pair we have a (A/Q) ' 16 correction, 
which is hardly distinguishable from a correction of order 1. 

Instead of truncating the perturbative expansion, we may achieve the same goal 
by putting a cut-off in the integral. In fact, consider the cut-off integral 

f ° dx i og 2* _L_ = / l0ST ^ dt t 2k e ~* . (3.15) 
Jo 1 — x Jo 

The saddle point of the integral is at t — 2k. If the saddle point is within the 
integration range, the integral is essentially the factorial of 2k, while for larger values 
of k it starts to grow like a simple power. Therefore, the cut-off acts like a truncation 
of the expansion. In order to have, as before, a truncation at k = l/(4a), we need to 
set the cut-off at log 1/(1 — x ) = 2k = l/(2a), corresponding to 

l-x = e~^. (3.16) 

In the worst case of production via gluon fusion in the MS scheme, we would have 
1 — xq = (A/Q) ' 32 . This leads to the rather paradoxical conclusion that the cut-off on 
the soft radiation should be imposed at values of Q much larger than A. For example, 
for the production of a 100 GeV object we would need a cut-off of the order of 14 
GeV. This would have to increase with energy. We observe that cut-offs of this kind 
are used in ref. J?|, in formulae (114) and (115). As a matter of fact, if we take the 
limit of large colour factor and small coupling of the formulae given there, we recover 
eq. (jjp . 

From the above derivation we see that the large corrections obtained have nothing 
to do with infrared renormalons, and it is easy to convince ourselves that they are 
a spurious effect. They were in fact not present in the original expression, eq. (|3.3| ). 
which is finite. It is also easy to show that the perturbative expansion of eq. (|3.3| ) 
has no factorially growing terms. This can be done in the following way. The k th 
coefficient of the expansion is given by the integral 

c k = — — % —±-\ g 2k N t~ n dN. (3.17) 

2TiiJc-ioo (N(N + 1)(N + 2)) 2 k\ & V ; 

We deform the integration contour as illustrated in fig. [I]. By choosing N sufficiently 



Figure 1: Integration contour for N for the determination of asymptotic 
properties. 



large, we get two contributions to our integral, from the circle and from the disconti- 
nuity along the negative axis. The integral on the circle is bounded by the expression 
C\og 2k Nq /k\ for some value of C. The discontinuity integral is instead given by 
(replacing iV — > — N) 



poo 

— / expfAHogr - 21og(iV(iV - 1)(N - 2))] Disc[log 2fe (-iV)] dN 

kl JNo 

36 r°° 

< — / exp [N log r + 2k log(7r + log N)] dN . (3.18) 

By saddle point integration of the right hand side we immediately see that the above 
expression cannot grow faster than 

1 /log2A;\ 2fc 



k\ y log I/7 



(3.19) 

Therefore the power expansion of eq. ( |3.3|) has an infinite radius of convergence, and, 
a fortiori, does not have factorially growing terms. 

Apparently, when performing the inverse Mellin transform to obtain the x-space 
expression of the cross section, we have simply thrown away subleading terms that 
would have compensated the factorial growth. In a leading-log sense we have every 
right to throw away subleading terms. However, if by doing so we generate unjustified 
factorially growing terms, we are certainly doing something wrong on physical ground, 
even if we are perfectly consistent with the leading-log approximation. 
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In the case of the full resummation formula, including the effects of the run- 
ning coupling, the above illustrated problem persists. We have in this case (see 
Appendix C) 

1 rC+ioo 

A(x) = — exp(\ogN gi (a s b \ogN))x- N dN 

Z7TI JC-ioo 

= —^(9(l-r]-x)exp[lg 1 (asb l))) x (1 + NLL terms) (3.20) 

in which an unintegrable singularity is met before x reaches 1, the Landau pole 
singularity. In the corresponding formula for the partonic cross section, neglecting 
NLL terms, we get 

£ dx exp[l gi (a s b l)}^C Q . (3.21) 

Expanding the above formula in powers of as, we would generate the same type of 
spurious factorial growth as found before. As before, the nature of the divergence is 
given by the integral 

/ dx exp[Z gi{otsbol)\ . (3.22) 
J o 

Using commonly available algebraic programs, it is easy to expand eq. ( 3.22|) up to 
large orders, and then study numerically the factorial growth. Expanding up to a s 2 
we have found the behaviour k\C^(boas) k , where C(k) is a slowly increasing function 
of k. If, for large k, C^) approaches a limiting value C, this corresponds to a power 
ambiguity of (A/Q) 2 ^ c . For gluon fusion in the MS scheme we get CV32) = 10.48, 
corresponding to a power ambiguity of (A/Q) ' 19 , while for qq annihilation in the MS 
scheme we get C^2) = 4.0, corresponding to (A/Q) ' 5 . These numbers are roughly 
consistent with those of the exact analysis performed for the fixed coupling case. 

We observe that, even if we modify the above x-space formula, by expanding it in 
powers of as and keeping only a fixed number of terms, the problem discussed earlier 
still persists. In fact, our discussion is relative to the case in which the exponent in 
formula (|3.20|) is expanded and truncated to order as- 

An x-space resummation procedure, similar to the one discussed here, has in- 
deed been adopted in the literature. In refs. [|, |(J a cut-off procedure is applied to 
screen the Landau pole singularity that manifests itself when x approaches 1. This 
cuts off both the divergence due to the Landau pole, and the spurious divergence we 
just described. We therefore argue that the uncertainties induced by this procedure 
are much larger than needed, since they introduce a divergence that is in fact not 
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present. In ref. 0, the Landau pole singularity is dealt with by using a principal 
value prescription. Subsequently, a cut-off is introduced in order to screen the large 
subleading effects that arise when performing the Mellin transform of the partonic 
cross section from N to x space. Our point is precisely that if these large subleading 
terms had been kept, they would have cancelled the factorially growing terms aris- 
ing from the leading terms after integration against the parton luminosities. In all 
these approaches, unphysically large cutoffs are needed in order to avoid the large 
corrections that arise at higher order in the perturbative expansion. 

We conclude that, when proposing a resummation formula for threshold effects, 
it is not enough to make sure that all leading terms are included in the formula. We 
must also make sure that we are not introducing subleading terms that grow very 
fast with the order of the perturbative expansion in the final physical result. In the 
next section we propose a resummation formula that is correct from the point of 
view of the threshold approximation, but does not induce any factorial growth in the 
perturbative expansion. 

4 The Minimal Prescription formula 

Our starting point is eq. ( j2.4|) . Its perturbative expansion has the form 

oo 

A N = J2 Ck(hgN)a k s (4.1) 

k=0 

where the coefficients Cfc(log N) are polynomials in log N. The resummed cross section 
can be formally written as a power expansion 

1 00 r C+ioc 

°r = E «I / f n(Q 2 ) c k (hgN) t~ n dN . (4.2) 

Observe that the integrals in the coefficients of the expansion are finite for C > 2, 
and no singularities occur at the right of the integration contour. The choice of C > 2 
is motivated by the usual Regge behaviour of structure functions, which implies that 
we cannot have any singularity in Fn(Q 2 ) to the right of the pomeron singularity, 
which is slightly above N — 1. The Landau pole for iV = exp 2a * 6o manifests itself in 
the fact that the series (|4.2j ) is not convergent. 

We now propose the following formula for the resummation of threshold effects in 
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the Drell-Yan cross section 

a res (r) = — / FUQ 2 ) &n(Q 2 ) t~ N dN, 2<C M p<N L = exp — — , 

(4.3) 

where A N (Q 2 ) is given in eq. (|2.8|) . The constant Cmp is chosen in such a way that 
all singularities in the integrand are to the left of the integration contour, except for 
the Landau singularities at N = Ni and N = iV|, which lie to the far right. We will 
call eq. (|4.3|) the "Minimal Prescription" (MP) in the following. Its justification relies 
on the following important properties, which will be proved in Appendix A: 

• The expansion ( f4.2| ) converges asymptotically to the MP formula. Observe that 
this would not happen if we had chosen a contour that passes to the right of 
the first Landau pole. 



The coefficients of the expansion (4.2) do not grow factorially. 



• If we truncate the expansion (4.2) at the order at which its terms are at a 

minimum, the difference between the truncated expansion and the full MP 
formula is suppressed by a factor 

e H a ? (4.4) 

where H is a slowly varying positive function. This suppression factor is stronger 
than any power suppression. 

We stress that with our MP formula we do not introduce any spurious factorial growth 
in the perturbative expansion. One may object that in this way we do not introduce 
any possible renormalon effect in the formula. Factorial growth due to renormalons 
is very likely to be present in the perturbative expansion. Our point is, however, that 
the leading terms in our expansion do not necessarily contain this factorial growth, 
and that renormalons present in a resummed expression therefore do not necessarily 
reflect the renormalons present in the full perturbative expansion. To be more specific, 
let us consider for a moment eq. (|2.4j ). It is clear that, if we perform the x integration 
exactly, we are indeed integrating over the Landau pole. However, since this formula 
is accurate at the NLL level at most, we may integrate it in the NLL approximation, 
and obtain formula ( |2.8|) , which has no trace of factorial growth. In particular, it 
was shown in ref. |14| that the leading IR arising from naively extending Eq. (|2.4j ) 
beyond the NLL level cancel in the full perturbative expansion. We take this result 



as a confirmation of the fact that the resummation of logarithmic effects at threshold 
does not teach us anything about the structure of power corrections. Resummation 
formulae should not, therefore, include any power correction. 

5 The Drell— Yan cross section 

We will not attempt to perform a detailed phenomenological analysis of resum- 
mation effects in Drell- Yan pair production in the present work. We will, however, 
assess the effect of resummation in the particularly simple example of the structure 
functions reported in Appendix B of ref. pi, for proton-antiproton collisions, and we 
will compare our results with the ones given there. In fig. 0we report the cross section 
as a function of r = Q 2 /S, normalized to the Born cross section. With this set of 
structure functions, which are Q 2 -independent, the .ff -factor depends only upon the 
ratio Q 2 /A 2 . The dashed curve is the NL cross section, including both the qq and qg 
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Figure 2: Drell-Yan pair production cross section in pp collisions, normal- 
ized to the Born result. 

subprocesses (given in ref. The dotted curve is the NL-resummed cross section, 
without the inclusion of the exact 0(as) result. The full curve is obtained by adding 
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to the NL cross section the NL-resummed contributions, after having subtracted the 
terms up to 0(as), that is to say 



g(t 



a 



L'». l nK T ) + T^T 



G 



+ 



Gy 



T - G r . 



T 



d 

Q S=0 c»q; 5 



<*S=0 



(5.1] 



Our result is consistent with the result of ref. |J, where analogous figures are given. 
There is however one important difference. The cut-off method used there to overcome 
the problem of the Landau singularity introduces an extra uncertainty, which is given 
by the spread of the various curves obtained using different cut-offs. This spread 
decreases roughly as 1/Q as Q 2 increases. Accounting for the fact that in ref. || also 
the so-called n 2 terms are exponentiated, our result is consistent with their band. 
It does not, however, agree with the central value, which is smaller in our case. 
Furthermore, even at the lowest energy, and for the very large value of r = 0.5, the 
fully resummed cross section is only 10% larger than the next-to-leading one. 



6 Heavy flavour production 



We will follow closely the notation of ref. ||21|| . The heavy flavour production cross 
section is given by the formula 

/•i / p H \ Am 2 

g= dx 1 dx 2 F(x 1 )F(x 2 ) g — — , p H = — tt , (6.1) 

J \X1X2J o 

where m is the mass of the heavy quark and S is the square of the total centre- 
of-mass energy. As before, for notational convenience, we have dropped here the 
parton indices. The scale dependence is also not shown explicitly. The partonic cross 
section depends also upon the heavy quark mass. Here we indicate explicitly only 
its dependence upon p H , which embodies its dependence upon the partonic centre- 
of-mass energy squared s = X\ x 2 S. In order to include the effects of soft radiation 
we parallel as closely as possible the approach followed in the Drell-Yan case. The 
leading logarithmic (LL) soft corrections in heavy flavour production have the same 



structure as in the Drell-Yan case [21, 0, since there are no collinear singularities 



arising from final state gluon radiation. We define 

— = [ dxi dx2F(xi)F(x2) 5(t — x\ x 2 )g ( — 
dr Jo V t 
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Ph\ 1 f c + io ° 



F 2 t~ n dN . 



t ) 27ci Jc- 

The inclusion of soft effects can now be performed as in the Drell-Yan case 



(6.2) 



da {res) 
dr 

Using now the identity 
and defining as usual 
we get immediately 



a 



n \ 1 rC+ioo 

pH ^ 1 f F 2 A^ F r~ N dN 



t J 2ni Jc-i 



f°° da 
a = dr — 
Jp H dr 



1 dz 



o z 



z N a(z) 



a 



(res) 



± [ C+i °° F* A» F *N-i dN. 

ZlXl JC-ioo 



(6.3) 

(6.4) 
(6.5) 

(6.6) 



The partonic cross section (after subtraction of collinear singularities) is 



&ij(s, m 2 ,p 2 ) 



fij(p,H /m ) 



(6.7) 



where p is the factorization scale (the renormalization scale is set equal to the fac- 
torization scale) and the dimensionless variable p is 



P 



Am 2 



(6.8) 



The functions fu have the following perturbative expansion 



Hp^ 2 /m 2 ) = ^\p)+g 2 {p 2 ) 



4/ ) (p) + ^ ) (p)ln^ 



(6.9) 



The lowest-order terms in Eq. ( |6.9| ) are explicitly given by ((3 = i/l — p 
f(o)/ _n 71 t rCf 



6 Nr. 



13 P (2 + p) 



(6.10) 



12 N2 - 1 



(4 + 4p - 2p 2 ) - In 



1 n 1+/3 



(3 1-/5 



4-4p 



+ C A 



3p 2 - In 



l+{3 



(3 1-13 



-4-5p 



(6.11) 
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and flj (p) = for all the other parton channels. 

The iV-moments of the expressions in Eqs. ( |6.10| ) and ( |6.11| ) follows 

vri T R C F r(AT+l) 



AO) 

Jqq,N 



AO) 

J gg,N 



3 

7T2 



N c T(N + 7/2) 
T R T(N + 1) 



(N + 2) 
1 



(6.12) 



x 



4 N 2 - 1 T(iV + 5/2) N + 3 
TV 3 + 9iV 2 + 20iV + 14 

2G P ; — : — 

(N+l)(N + 2) 



C A 



N 2 + 8N + 11 
2iV + 5 



(6.13) 



The resummation effects are embodied in the factor A^r. While in the Drell-Yan 
case this resummation factor is only associated with the qq subprocess, in the heavy 
flavour case both qq and gg subprocesses are involved. The explicit expression of A^v 
will therefore depend upon the subprocess. To leading logarithmic accuracy we have 
(a s = a s (m 2 )) 



lnA^(m 2 



In Ng ijA (b a s In N) + 0(a k s \n k N) , 



(6.14) 



where the functions g^ \ are related to the function gi in Eq. (|2.10|) by simple colour 
factors. More precisely we have 



9gg,lW = 7T 9iW ■ 



(6.15) 



The factorization scheme dependence is of course contained in g\. We will therefore 
take eq. ( |6.6D with C = Cmp as our MP formula for heavy flavour resummed cross 
sections. 



7 Heavy Flavour cross section: phenomenological results 

In this section we present some phenomenological applications of the resummation 
formulae presented above. To start with, we present in fig. [5] the partonic cross 
sections for production of a pair of heavy quarks of mass rriq = 175 GeV, plotted as a 
function of rj = (1 — p)/ p. The figures show the Born, NLO and resummed results for 
both the gg and qq initial states. The resummed partonic cross section is defined as 
the sum of the contributions of order a| and higher from eq. (|6.6|) and the fixed order 
NLO result. These figures can be compared with similar ones in ref. 0. Notice that 
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while in that work the growth of the resummed cross section at small rj is damped by 
a cut-off, in our approach the growth at small rj is automatically controlled. Notice 
also that the resummed cross section is a smooth function of p down to values of 1 — p 
of the order of 10 -3 , which is of the order of the ratio A/tuq. Since the energy of the 
soft radiation is of the order of (1 — p)mQ, at these values of p it is numerically of 
the order of A. This is the point at which the Landau pole is expected to influence 
the results, and non-perturbative physics to set on. We interpret this behaviour as 
a confirmation of the correctness of our procedure. It would make no sense to cut 
off the partonic cross section at values of p corresponding to soft gluon radiation of 
several GeV. For 1 — p < 10" 3 the resummed cross section starts to oscillate, but it 
remains integrable. 

In all our phenomenological studies of resummation effects in heavy flavour pro- 
duction we have used the structure function set CTEQ1M ||22|| . The importance of 



the resummation effects is illustrated in figs. f|, | and |, where we plot the quantities 

^gg ^qq ^gg + ^qq 

r (gg) ' ' Agg) , Am) 



(7-1) 

^NLO "NLO ^NLO + °NLb 

Here S is equal to our MP resummed hadronic cross section in which the terms of 
order a 2 s and a| have been subtracted, and <7(nlo) is the full hadronic NLO cross 
section. The results for b at the Tevatron can be easily inferred from fig. ^, since the 
qq component is negligible at Tevatron energies. 

For top production, we see that in most configurations of practical interest, the 
contribution of resummation is very small, being of the order of 1% at the Tevatron. A 
complete review of top quark production at the Tevatron, based upon these findings, 
has already been given in ref. |23fl . We also observe that, for top production at the 



LHC, soft gluon resummation effects are negligible. Of course, in this last case, there 
are other corrections, not included here, that may need to be considered. Typically, 
since the values of x involved are small in this configuration, one may have to worry 
about the resummation of small-x logarithmic effects [J24| . 

We see from the figures that in most experimental configurations of interest these 
effects are fully negligible. One noticeable exception is b production at HERAb, 
at v^S* = 39.2, where we find a 12% increase in the cross section. This correction is 
however well below the uncertainty due to higher order radiative effects. For example, 
from the NLO calculation with the MRS A' |25| parton densities and = 4.75 GeV, 



we get a b i = 10.45^'^ nb, a range obtained by varying the renormalization and 
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Figure 3: Partonic cross section for the production of a 175 GeV heavy 
quark pair, qq initial state (left) and gg initial state (right). The dotted lines 
are the Born result, the dashed lines the NLO result and the solid line the 
resummed result. Af s =152 MeV. 
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Figure 4: Contribution of gluon resummation at order a| and higher, relative 
to the NLO result, for the individual channels and for the total, for bottom 
production as a function of the CM energy in pp collisions. 
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Figure 5: Contribution of gluon resummation at order ag and higher, rela- 
tive to the NLO result, for the individual subprocesses and for the total, as a 
function of the top mass in pp collisions at 1.8 TeV. 




Figure 6: Contribution of gluon resummation at order Og and higher, relative 
to the NLO result, for the individual channels and for the total, as a function 
of the CM energy in pp collisions. 
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Figure 7: Contribution of gluon resummation at order a s and higher, relative 
to the NLO result, for the individual channels and for the total, for charm 
production as a function of the CM energy in pp collisions. 



factorization scales from m&/2 to 2m&. Thus the upper band is 80% higher than the 
central value, to be compared with a 10% increase from the resummation effects. This 
result is much less dramatic than the results of ref. IF2D . 



In fig. [?| we also show the effect of resummation in charm production. Due to 
the large uncertainties that plague charm production |p6| , this plot should only be 
considered for orientation. 

We conclude with a few remarks and checks about our result. We will focus on 
top production by qq annihilation, with m t = 175 GeV at the Tevatron. First of all, 
we checked that the full resummation formula that we use is well approximated by 
its expansion in powers of as- We define 

M 

4 es) = E°i. ( 7 - 2 ) 

k=2 

where crjt is the contribution of order a k s . We computed each term in the expansion 
up to order a%. The results are displayed in table [I]. We see that the convergence 
properties of the expansion are extremely good, and, up to the order we have probed, 
there is no sign that we are near the breakdown of the expansion. As a second 
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observation, we notice that the term 03 in our resummed formula is about 7% of the 
Born result. The full 0(a\) correction is instead 20% of the Born term. It is easy to 
illustrate the contribution of the 03 term to the partonic cross section. We start with 
the full NLO cross section written in the form of eqs. ( |6.7| . |6.9| ) 

{NLO \p) = ^\^{p) + ^a s h qq { P )) (7.3) 



a 



where 



77T 



h qq { P )=fV{ P ) + f^{p)\og^- 2 . 



(7.4) 



The functions f q q , f q ^ and fj^ are defined in ref. and p is the factorization and 



renormalization scale. It is easy to show that the truncated resummed result at order 
0(a s ) can be obtained by using the following partonic cross section 



4 es \p) = ^ 2 (fS\p) + ^sK, 

fit 



(7.5) 



where the function h qq can be obtained by expanding the resummation function 
obtains 



A qq F N (m 2 ), as given in eq. (|6.14j ), up to order as- After some simple algebra one 



KM 



c F 

7T- ./,, 



■ 2 1 dyfS\p/y) 



logl/y 

The meaning of the plus distribution is as usual 



(log log l/y + 7 S ) 



(7.6) 



dy[G(y)} + F(y)= / dy G(y) (F(y) - F(l)) . (7.7) 
Jo 

Notice the presence of the subleading term proportional to je in the equation. This 
term is cancelled if we use, when performing the x integration in eq. Q2.4p , the relation 
1 — x N = 9(1 — x — ^-jf-) (which is accurate up to NLL |J) instead of the approximate 
9(1 — x — -k). This amounts to the substitution: 



1 N 

one, 1 — x 



In A^ N (m 2 



1 + 7 *5lulv) ln 



(7- 



Equation (|7.6|) then becomes 



KM 



c_ F 

7T- 



2 , dyfgf(p/y) 



log l/y 



log logl/?/ 



(7.9) 



In fig. H we give the contribution of the h q q (p) compared with the quantity that 
corresponds to the full NLO correction. 
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3.71 


0.256 


8.27 x 10~ 3 


-3.24 x 10~ 4 


-5.88 x 10~ 5 




0.066 


2.0 x 10~ 3 


-9.63 x 10~ 5 


-1.48 x 10~ 5 


4.8 x 10" 8 



Table 1: Top pair production at the Tevatron via the qq channel, for mt = 
175 GeV and CTEQ1M parton densities. Order-by- order contributions to the 
fully resummed formula (first line), and accuracy of the truncated perturbative 
expansion relative to the fully resummed result (second line). 
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Figure 8: Comparison of the partonic cross section at 0{o? s ) for the MP 
resummation formula (solid line), for the resummation formula with the term 
proportional to 7^ included (dotted line), and for the exact 0(a|) formula of 
ref. with [i = m t /2, m t and 2m t (dashed lines). 
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We plot the full NLO correction for \x = mt/2, mt and 2mt. We see that the 
h' q q term is consistent with the exact next-to- leading result, given the spread of the 
latter induced by the renormalization scale dependence. According to the choice of 
subleading terms, the 03 term can substantially change, and can be brought to almost 
coincide with the exact result for \x = m t . In our case, for example, inclusion of the 
term proportional to 7^ in formula ( 2.16|) would cancel exactly the corresponding 
term appearing in formula (|7.6| ), which then reduces to Eq. ( |7.9|) . This last formula 
gives a result that is very close to the exact one for = mt in the most important 
kinematic region of p ~ 1, as can be seen from fig. ||. This is in agreement with the 
findings of ref. |§, where it was pointed out that the LL truncation of the 0(a%) 
terms provide a very good approximation to the full NLO cross section, as long as 
fi = m t . However it should be stressed once more that this result is accidental, as it 
would not hold for a different choice of renormalization scale. 

Notice that altering the structure of the subleading terms in the exponent of the 
Mellin-space coefficient function will not affect the asymptotic properties of its per- 
turbative expansion and the integrability of its a;-space MP transform. We explored 
the numerical impact on the contributions of order a| and higher of including the 
subleading terms proportional to 7^ in the exponent of the coefficient function. In 
the notation of Eq. we get 5 qq -/ cr NL = 0.013, an effect which is larger than the one 
found previously, but still negligible. The difference from the result of fig. |5| should 
be taken as an estimate of the uncertainty coming from the unknown next-to-leading 
logarithmic terms in the exponentiated coefficient function. As such, it is a purely 
perturbative uncertainty, which cannot be separated from the uncertainty due to the 
change in renormalization scale or factorization scheme. The size of these latter un- 
certainties, estimated to be of the order of 10% PB[ , can consistently accommodate 
for the 1% effect we found. 



8 Jet Cross Sections 

In this section we present, as an additional application of our formula for the soft 
gluon resummation, a study of corrections to the invariant mass distribution of jet 
pairs produced in pp collisions at 1.8 TeV. The interest in the effects of resummation 
on the behaviour of jet cross sections at large energy is prompted by the discrepancy 
between the single-inclusive jet-px distribution at large pp, as measured by CDF [^7j 
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and the result of the NLO QCD predictions |p8| . For simplicity we will study the 



effects of soft gluon resummation on the invariant mass distribution of the jet pair, 
which is, from a theoretical point of view, very close to the Drell-Yan pair production. 
Observe that other distributions, such as the pt of the jet, have a rather different 
structure from the point of view of soft gluon resummation. In fact, while the jet pair 
mass is only affected by the energy degradation due to initial state radiation, the pr 
of the jet may also be affected by the transverse momentum generated by initial state 
radiation, and by the broadening of the jet due to final state radiation. 

A study of the jet pair mass distribution is not of purely academic interest, since 
also for this variable an analogous discrepancy between data and theory has been 
observed ||29|| . Studies of resummation effects in the inclusive spectrum of jets are 
in progress (M. Greco and P. Chiappetta, private communication). 

Before presenting the results, we briefly discuss the key elements of the calcu- 
lation. Contrary to the case of Drell-Yan and heavy flavour production, light jets 
are produced at Born level from all possible initial states, gg, qg, qq and qq. The 
resummed invariant mass distribution of dijets, at LL order, is therefore given by 

~^-= E F^f** A£. iJV t~ n dN , (8.1) 

CLT . ._, _ -i Z711 JC—ioD 



where s = St = m 2 jj. The functions A^- )JV are given by eqs. (|6.14j ) and (|6.15|) , 
supplemented by 



lnA^^mi,) = \nN g ggtl (b a s \nN) + O(a k s \n k N) (8.2) 

0«,i(A) = %^i(A). (8.3) 

The partonic Born cross section &^ (s, Q* min ) depends on the range of integration for 
the partonic centre-of-mass scattering angle 9* (6* is related to the rapidity difference 
V* = (Vi ~~ ? ?2)/2 of the two jets by sin#* = 1/ cosh?]*). To avoid the Rutherford 
singularity, we will keep this angle strictly larger than zero. While the absolute 
production rate depends on the choice of Q* min , we will now show that the i^-factor, 
i.e. the ratio of the resummed cross section to the LO one is to good approximation 
independent of it. Therefore the results for the ii'-factor will be rather independent 
of the details of the experimental cuts. It is a well-known fact [3(| that the LO 



amplitudes for parton-parton scattering processes are related to one another, in the 
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small scattering angle limit, as follows 

4 /4\ 2 

da gg : da qg : da qq = 1 : - : , (8.4) 

where the indices refer to the pair of partons in the initial state, q being an arbitrary 
quark (or antiquark) flavour. Even at 90° this approximation is good to about 10%, 
becoming better and better when the cross sections are integrated over a larger and 



larger range of the scattering angle. Using eq. ( |8.4j ), we can rewrite the expression 
for the resummed jet invariant mass cross section as follows 

^(res) i r C+ioo 

MMiUTr-/. H 2 N r~ N dN, (8.5) 



dr yyy 1 mm ' 2m Jc 
with 

H N = F 9tN N^ 9l{X) + 1 £ F i>N N^ . (8.6) 

The .fT-factor is therefore given by the following expression, in which all dependence 
upon the acceptance cut 9* has disappeared 



J^-(res) 1 rC+ioo 

' ' Hlr- N dN 



dr _ 2m J c x ' 



j^-(LO) 1 r C+ioo 

da 1 f E 2 N r~ N dN 



5.7) 



dr 2m Jc- 

E N is the Mellin transform of the standard effective structure function, defined by 

4 

En = F 9) n + — F^n . (8.8) 

While a more accurate implementation of the experimental acceptance will eventually 
be necessary^, our approximation is more than adequate to provide a first estimate 
of the resummation effects. 

In fig. we show the following quantities 



<j(3) §(3) x(3) 5 (3) , § (3) , 

U gg "qg °qq "gg ^ "qg ^ "qq 

a (2) ' a (2) ' a (2) ' ^(2) 



.9) 



where 5^ is equal to our MP resummed hadronic cross section in which the terms of 
order a 2 s have been subtracted, and = aj® + + is the full hadronic LO 
cross section (of order a|). We use as a reference renormalization and factorization 
scale for our results fi = Mjj/2. Notice that for large invariant masses the effects of 
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Figure 9: Contribution of gluon resummation at order a|. and higher, rel- 
ative to the LO result, for the invariant mass distribution of jet pairs at the 
Tevatron. 




Figure 10: Contribution of gluon resummation at order ctg and higher, 
relative to the truncated 0{a? s ) result, for the invariant mass distribution of 
jet pairs at the Tevatron. 
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higher orders are large. To understand how much is due to the first order corrections 
(which are exactly calculable [|3T]) and how much is due to corrections of order a\ 
and higher, we show the following quantities in fig. |K] 

tf(4) g(4) x(4) ^(4) , 6 (4) , 5 (4) 

"gg "qg °qq "gg ^ "gg J_ "qq f8 in"! 

a (3) ' a (3) ' a (3) ' a (3) ' ^°- ±u ^ 

where 5^ is now equal to the MP resummed hadronic cross section with terms of 
order a| subtracted, and is an approximation to the full NL cross section, summed 
over all subprocesses, obtained by truncating the resummation formula at order a|. 
This figure shows that indeed most of the large K factor is due to the pure NLO 
corrections, with the resummation of higher-order soft gluon effects contributing only 
an additional 10% at dijet masses of the order of 1 TeV. Removing the subleading 
term proportional to Je from the exponent of the coefficient function, in the spirit of 
the discussion at the end of the previous section, we found only a slight increase of 
the contribution from the terms of order a| and higher (an additional absolute 5% 
for dijet masses of 1 TeV). 

As we will discuss in more detail in the next section, these results should only be 
taken as an indication of the order of magnitude of the correction, since we have not 
included here a study of the resummation effects on the determination of the parton 
densities. From this preliminary study it seems however unlikely that the full 30-50% 
excess reported by CDF for jet pt's in the range 300-450 GeV could be explained by 
resummation effects in the hard process. It is possible that the remaining excess is 
due to the poor knowledge of the gluon parton densities at large x, an idea pursued 
by the CTEQ group [|3"2]| . This problem has been studied also in ref. [Q. 

9 A few additional remarks 

A fully consistent treatment of the resummation effects requires the use of parton 
densities that 

(i) are extracted from low-energy data by taking into account the resummation 
effects for the corresponding scattering process; 

(ii) are evolved in Q 2 using resummed anomalous dimensions. 

4 For instance, at small 9* a more refined resummation procedure would lead to the replacement 
^jj ~ * ^jj sm2 ($*/2) in the hard scale of the resummed coefficient functions - N in eqs. 
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In both steps (i) and (ii), the Minimal Prescription should be implemented. 

In the case, for example, of jet production, it is quite possible that resummation 
effects significantly influence the determination of the large- a; structure functions from 
low energy data. The Tevatron jets of the highest energy probe the partonic densities 
at values of x of the order of 0.5, which one could argue are far enough from the x — > 1 
region to make the Sudakov effects small. However, the large values of Q 2 involved (of 
the order of 10 6 GeV 2 ) are such that the Q 2 evolution of the parton densities from low 
energy is significant, and we are therefore sensitive to the input structure functions 
measured at low Q 2 in regions of x significantly larger than 0.5. It would therefore be 
important to reexamine the extraction of the large- x non-singlet structure functions, 
in the light of the resummation results for the DIS process, before firmer conclusions 
can be drawn on the significance of the present jet cross section discrepancy. 

In the Q 2 evolution of the parton densities, one has to use the equation 

= E^wMQ 2 )) F hN {Q 2 ) , (9.i) 

with anomalous dimensions jij,N that include the resummation of logarithmic contri- 
butions to the same accuracy as in the coefficient function. For instance, the analogue 
of Eq. (|2.8|) is the following expansion 

TWv(as) = 7§ } (&o« 5 In N) + a s if {b a s In iV) + . . . . (9.2) 

The leading and next-to-leading resummed anomalous dimensions 7}* and 7-? 
have, in general, singularities related to the Landau pole and these singularities (like 
those in Eqs. ( f2.10| ), ( |2.11| ), and ( [2.15| ), Q2.16D ) depend on the factorization scheme 



in which the parton densities are defined. Thus, for the purpose of implementing the 
Minimal Prescription, it is particularly convenient to use the MS definition of the 
parton densities. Indeed, in the MS factorization scheme the resummed anomalous 
dimensions have no singularity associated with the Landau pole and, more precisely, 
they have the following simple form H, EM 



7.SvK) = ~ A(a s ) + 0(a%)\\nN + 0(l) , (9.4) 
where A(as) is given in Eq. (p~5|) and j]f s N (as) = 0{1/N) for N — ► 00 if i ^ j. 



^ S n(*s) = ^ A(a s ) + 0(a s ) lnN + 0(l) , (9.3) 

,MS 
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Note, however, that the main feature of the Minimal Prescription, i.e. the ab- 
sence of factorially-growing coefficients, remains valid in any factorization scheme. 
The only difference in using different schemes is in the fact that one can introduce 
non-perturbative corrections which, although suppressed by more than any power of 
A/Q at fixed 1 — r, can actually have a different overall magnitude. As long as Q 2 
is sufficiently perturbative and r sufficiently far from the hadronic threshold, this 
difference should not sizeably affect the actual value of the cross section. Obviously, 
approaching the essentially non-perturbative regime (Q 2 — ► 1 GeV 2 ,r — ► 1), a phys- 
ically motivated treatment of non-perturbative effects has to be introduced. This is 
beyond the scope of the present paper. 



10 Conclusions 

In our study of resummation procedures for partonic-threshold corrections, we 
have found that the factorial growth due to the infrared renormalon is only a minor 
problem, when compared with the large spurious factorial growth that is generated 
when the x space resummation formulae are limited to the LL or NLL level. Although 
we consider the recent indication of the cancellation of 1/Q renormalon effects in 
Drell-Yan pair production [111 as an important progress, we find this spurious factorial 



growth to be a much more serious problem in practice. In fact, while a 1/Q power 
ambiguity in the Tevatron top production cross section should be below the per cent 
level, the ambiguity induced by the spurious factorial growth is at the level of 10%. 

In the present paper we have proposed the Minimal Prescription for the resum- 
mation of partonic-threshold effects in hadronic processes. This formula has a per- 
turbative expansion free of factorially growing terms. The ambiguity arising from 
its asymptotic nature is in fact exponentially suppressed, behaving as exp(— H(l — 
x)Q/A). We would also like to remark that certain kinematic constraints are respected 
by the MP formula. For example, for N = 1 the MP formula gives no resummation 
corrections. In general, sum rules associated to low moments receive small corrections 
in our procedure. 

As far as our phenomenological results are concerned, we have found that in heavy 
flavour production, in current experimental configurations, resummation effects are 
negligible. The process of high mass dijet production at the Tevatron, being much 
closer to threshold, has instead non-negligible corrections. We wish to remind the 
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reader that, in this last case, in order to perform a reliable phenomenological predic- 
tion, resummation formulae should be applied not only to the production process in 
question, but also to the processes that have been used to extract the parton densities 
and to their evolution equations. 
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APPENDIX A: Asymptotic nature of the resummation for- 
mula 

In this appendix we will study the asymptotic nature of our resummation formula. 
We are dealing with the LL function 

1 rC+ioo 

G(a s ) = - — : / dN r~ N h N (X) , h N (X) = expflog AT/(A)] , A = a s b logTV. 

(A.l) 

Observe that we will consider h N (\) as an independent function of N and A. We will 
focus upon the MS case, in which 

/(A) = j [(1 - 2A) log(l - 2A) + 2A] (A.2) 

with c = / (tt&o), but the proof can be easily generalized, since it only relies upon 
the analyticity properties of /(A), and the fact that it does not grow too strongly at 
infinity. Because of the same reason, the proof can also be generalised to subleading 
logarithmic accuracy. 

We introduce the notation 

n \i 

h N (\)\ n = £/$(o)-. 

i=o l - 

G(a s )\ n = J dN t~ n h N (\)\ n ■ (A.3) 

It is clear that G(as)\ n is the truncated expansion of G(as) to order n in as- We 
will show that G(as)\ n is asymptotic to G(as), that is to say that 



SG(a s ) = G(a s ) - G(a 5 )|n-i = 0(a n s ). 



(A.4) 
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We begin by stating a few useful results for analytic functions. For a generic analytic 
function F(z) we have the identity 



n-l 



dxi 



x 1 



dxo ■ ■ ■ 



dx n F^ ^(x n ) 



(n-l)! Jo 



,z — x 



\ n - l F {n) (x) dx. 



(A.5) 



where we assume that the integrations are performed along a straight path from to 
z. We then use Cauchy's inequality, which states that for a function F(z), analytic 
in a circle of radius r around the point z we have 



|FW(* )|<^max |F(*)| ■ 

T n \z-z \<r 



(A.6) 



Assuming that F(z) is analytic in a circle of radius R around z = 0, from eq. ( A.6|) 
it follows that 



\ Fin \z)\ < ,„ 'i I \ max \F(w) 

' W ' _ (R- \z\) n \w\<R l V 1 

which together with eq. ( |A.5| ) yields the following bound 



n-l 



F(z)-J2FM(0) 



k=0 



< n max \F(w) 

\w\<R' 



2 I (\z\ - x) 



n-l 



(fl-x) 



dx 



(A.7) 



(A.8) 



which holds for \z\ < R. We now use the relation 



(t-y) 



n— 1 



o (i?-Z/) n 



M n f 1 (i-y) 



n-l 



o (1 - yt/.R) n 



< - ( 4 



i? 



(A.9) 



(where the last inequality becomes an identity in the limit n — > oo) and get the bound 



n-l 



F(z)- 5> w (o)tt 

fc=0 



< max 



i? 



|to|<ii 



(A.10) 



which we will apply in the following. 
Consider now our expression for SG 



SG(a s ) 



2ni Jc 



C+ioo 



dNr~ N (MA)-MA)| n _!) 



(A.ll) 



We will apply the bound ( A. 10 ) to h^(X), considered as a function of A at fixed 
N. This function is analytic up to the Landau pole at A = 1/2, so that in our case 
R — 1/2. We bent the N contour towards the negative real axis, as depicted in fig. [TJ 
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with Nq of order 1. Choosing any real < r < 1 (close to 1) we split our integral 
into a near (| logiV| < r/(2as&o)) and far (| logiV| > r/(2as&o)) part, 



G{a s ) = G K (as) + G > {a s ) 
G<>(as) = 77— / 

2lTt J\logN\<>r/(2a s b ) 

For the near part, using eq. ( |A.10|) , we obtain immediately 



dN r~ N h N (X) 



(A.12) 



i<jg< 



|logAT|<r/(2a s 6o) 



dNr~ N (h N (X)-h N (X)\ n ^) 



< 



where 



and 



|A r |<e r ''( 2Q! S i, o) 



9n 



t _ n exp[C|logiV|] 



\dN r- N 



1 — r 
exp[C|logN|] 



1 — r 



(2a s b \log N\) n <g n a n s , (A.13) 
(26 |logiV|) n (A.14) 



C= max |/(A)|. 

|A|<r/2 



(A.15) 



By inspection it is easy to see that the far part is exponentially suppressed. In fact, 
we can write 

|<JG>| = / dN r- N (h N {\) - MA)L-i) 

J|logAf|>r/(2a s fe ) 

< / \dN\e- lo ^W \h N (X) - h N (X)\ n ^\ . (A.16) 

J\logN\>r/(2a s b ) 

Since f(X) is logarithmically bounded in the far region, h^(X) grows at most as 
(N) lnlnN . Instead, h^(X)\ n is polynomial in log AT. Therefore, the dominant factor in 
the integrand is the exponential, which is of the order of exp(— exp(r/(2as&o)) log 1/r). 
We have therefore proved that the MP formula is asymptotic to its formal as expan- 
sion. 

It is easy to find, using saddle point methods, the rate of growth of g n . We find, 
for n <C log 1/r, the following power behaviour 

c V 



g n « 2b log ■ 

V log 1/r y 

and for n 3> log 1/r we obtain the coefficients 

n 



g n « 2b log 



log 1/r 



(A.17) 



(A.18) 
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whose growth is faster than any power, but much slower than factorial. We can 
estimate the value of n at which a$g n has a minimum from the equation 



a n s +1 9n+i 



log 



n+l 
log 1/r 



2a s b log 



n + l 
log 1/r 



log . 



a 's9n V l0g ^iV 

The value at the minimum is obtained from 



^77 2asbo log 



n + l 
log 1/r 



a n s g n « exp \n\og 
by using eq. ( |A.19[ ). We get 



2a s b log 



n 



log l/i 



n 



a s g n « exp 



log 



log 1/r 



The value of n at the minimum is 



n 



e 2a s b o log 



T 



which inserted in eq. ( A.21 ) yields 

a n s g n « exp ( -2a 5 6 e~s° 



log- 



= 1 . 
(A.19) 

(A.20) 

(A.21) 

(A.22) 
(A.23) 



Remembering that as bo = 1/ log(Q 2 /A 2 ) we see that, taking r — > 1, these suppression 
factors correspond to exponentials of the form (neglecting logarithmic factors in the 
exponent) 



-H 



QQ-r) 



e " a , (A.24) 

which is suppressed by more than any power of Q. Notice that as Q(l — r) approaches 
A the above expression becomes of order 1. In fact, in this limit, the total mass of 
the radiation that accompanies the production of the object with mass Q becomes of 
order A. This is clearly a regime over which we do not have any perturbative control. 

As a last point, we discuss the inclusion of other A^-dependent factors in the 
integrand like the partonic cross section in the case of heavy flavour production, or 
the parton luminosities. Partonic cross sections have a negative-power behaviour in 
N for large N, so they will appear in the integrand as factors of exp[— MlogA], 
which simply modifies the value of the constant C. Similar effects would be given 
by common parametrizations of structure functions. Of course, we cannot guarantee 
that the structure functions themselves do not have a Sudakov like behaviour at large 
x. The only statement we can make is that in the MS scheme, if such a behaviour is 
not present in the initial conditions, it will not arise because of evolution (see sec. 0). 
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APPENDIX B: Numerical implementation of the MP for- 
mula 

In principle, there are no difficulties in the numerical implementation of the MP 
formula. However, the modern sets of parton densities are usually provided in terms 
of numerical codes or of parametrized expressions which, in practice, cannot be used 
to evaluate their iV-moments in analytic form for arbitrary complex values of N. To 
overcome this practical difficulty we have to rewrite the MP formula in terms of an 
x-space convolution of the parton densities and the inverse Mellin transformation of 
the coefficient factor Ajy. Since A N (Q 2 ) in formula Q has singularities to the right 
of the integration contour, much care has to be taken when turning formula (|4.3Q 
to an a;-space integral. We begin by observing that, for large N, A^(Q 2 ) (in the 
MS scheme) is suppressed by more than any power of N as N —>■ oo. We will also 
assume that Fn is suppressed by some powers of N as iV — > oo. We will limit the 
present discussion to the Drell-Yan case, although the extension to other cases is 
straightforward. 

We rewrite formula fl4.3| ) by setting 

rl dx 



where 



m<r)= / —x £&Q) (b.i) 

JO X 

£(*,Q 2 ) = / F(x/z)F(z)-. (B.2) 

J x Z 



We obtain 



O T 



1 ,c MP +ioo r ,1 dx 

— / t A N (Q ) / — x L{x,Q 

2m Jcmp-ioo [Jo x 

1 dx „, „n N 1 rCMP+i°o / T \- N 



dN 



£(x,Q' 2 ) 7T-. I' I") &N(Q 2 )dN. (B.3) 

o x 2m Jc M p-ioo \xj 

Observe that the x integration extends from to 1, not from r to 1. This is because 
of the fact that, due to the Landau singularity to the right of the contour, the Mellin 
transform of A^(Q 2 ) does not vanish when its argument is greater than 1. It does 
however vanish very fast when its argument is above 1 by an amount greater than 
A/Q, so that the basic parton model assumptions are not violated. Let us in fact call 
z = t/x. For z > 1 we have 

— / A N (Q 2 ) z~ N dN =- e- m °^ lmA N+le (Q 2 ) dN . (B.4) 

2m Jc M p~ioo 7T JN L 
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Using the formula 



-Nt 



E 

fc=0 



(B.5) 



where 5^ h \t) is the k th derivative of the 5 function, we obtain 



2m JC] 



Cmp+*°° 



A N (Q' 2 ) z~ n dN 



M p -ICO 



5{\ogz)- 

71 JN L 



T A /^2N dN 

lmA N+ie (Q 2 ) — 



dN 



o&z)~ I lmA N+ie (Q 2 ) + • 

71 JN L N Z 



(B.6) 



Using saddle point arguments it is possible to show that the first term is of order 1, 
while the subsequent terms are suppressed by powers of A/Q. Let us in fact consider 
the integral 



1 

71 JN L 



lmA N+ie (Q 2 ) 



dN 1 fCmp+ioo „ dN , 

mL-W^JF*- (B7) 



The integration can be turned into a A = OJ560 log N integration, where the contour 
C comes from 00 — ei, encircles clockwise the Landau singularity at A = 1/2, and then 
goes to 00 + ei. In the MS scheme we have 

1 1 



h 1 • / ex P V u 
o as 2,7x1 Jc \olsO[ 



/ r 



m 



2A A + (1-2A) log(l-2A) 



dX 



(B.8) 



where r 



/ ipQTx) . The integration contour can be deformed into a steepest 

1 1 1 ~ ~ 1 and the saddle 



descent contour, where the saddle point is located at A = ^(1 — e 2 
point estimate of the integral is 



exp 



e 4r 



otsh 



1 - 



m 
2r 

exp 



e 2r 



1 



has 2m 
m 



exp 



( 2r m - 2 



\a s b 



e~r X z dX 



2r 



— e 2r 



(B.9) 



We have verified numerically that the above estimate is in fact quite good, although 
terms down by a power of as have been neglected in the saddle point approximation. 
It is apparent from the formula that, except for m = 0, I m is power-suppressed, with 
an exponent increasing with m. In our philosophy, however, all the terms of the 
expansion should be kept, since we have chosen not to include any power- suppressed 
effects in our definition. 

In principle there is nothing wrong in the integral ( p.3|) , it is convergent, and 
it should be possible to perform it numerically. In practice, however, this is not 
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convenient. Even if a careful importance sampling is implemented in order to correctly 
treat the region around z — 1, the integrand performs a few oscillations, in which 
large cancellations occur, so that a direct numerical integration proves in this case to 
be too time-consuming. In order to overcome this problem we use the following trick. 
We introduce a fake luminosity C^°'(x,t), defined as 

£(°\x, r) = A (1 - x) a x- 13 (1 - P{x)) . (B.10) 

where P(x) is a polynomial in x. We then define 

Xl =r + (z- 1)77, U = C<®(xi,Cf), (B.ll) 

where i — 1, ... ,4 and 77 is a small quantity, usually taken to be of the order of A/Q. 
We then choose 

P(x) = B (x — X\) (x — x 2 ) (x — £3) , (B.12) 
and A, B, a, (3 in such a way that 

C®(xi,T) = li, z = l,...,4. (B.13) 

These conditions can be easily solved to yield 

log - log log — log — 

h x 2 h x 3 



a 



P 



A 



dog 


Xi 




X2 


1 - 


Xi 


1 - 


£3 


1 - 


Xi 


1 - 


x 2 


1 - 


Xi 



Xi 1-Xi 1-xi Xi 

log — log log log — 

x 2 1 — x 3 1 — x 2 X3 



Zi 1-xi Zi 1-xi 
log - log log — log 

h 1 ~ x 2 h 1 ~ 3^3 

Xi 1 - Xi 1 - Xi Xi 

log — log log log — 

x 2 1 - x 3 1 - x 2 x 3 

Zi 



;i - xi) q x/ 



5 = — ^ 4; 7, r . (B.14) 



A(l-x 4 ) a x± p - I 
A (1 — x 4 ) Q x~l^ (x 4 — Xi) (x 4 — x 2 ) (x 4 — x 3 

We now rewrite formula ( |B.3| ) in the following way 

r 1 dx f „, ,\ 1 /-CMP+ioo f T \~ N 

cr(r) 



JO X V 7 Z7TZ JC M p-*oo Vx/ 

+ t— / r"" C^(r) A n (Q 2 ) dN , (B.15) 

Z7TI jCmp^*oo 
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where as usual 

£ ( n\t)= f -x n £(°\x,t), (B.16) 
Jo X 

which is easily performed analytically for possibly complex values of N. The difficult 



region of x — > r is now strongly suppressed in the first integral of eq. ( |B.15| ), and the 
whole formula can be implemented numerically with no further difficulties. The only 
problem one encounters is instabilities with structure function packages that perform 
piecewise interpolation of tables, since they do not give a smooth structure function. 
Thus, we found it much easier to use the CTEQ sets, which are analytic in x, rather 
than the MRS sets, which use an interpolation procedure. 



APPENDIX C: Some formulae and theorems on the Mellin 

transforms 

Let us begin by discussing the basic result 

/ ^-iVlogz+logiVFCaslogiV) = ( 1 _ x ) e lF(a s t) x [l + E(a S , I)) (C.l) 

where 

I = — 1ok(— losrx) ~ lo£ 



1 — X 

F{z) = Y.hz k 



oo 

k 



k=l 

oo k 

E(a s ,l) = X>!E e ^'- (C2) 

k=l j=0 

Therefore, the E term is a NLL correction, since its power expansion has no terms 
with more powers of I than of as- Equation (|C.1| ) should be interpreted in the sense 
of its formal power expansion in as- A simple proof has been given in ref. JOJ. For 
completeness, we shortly report here the basic argument. We rewrite the integral in 
terms of z — — iVlogx 

e -N\o Z x+\o % N F{a s \o^N) _ f g2+[log z+l] F(a s [log z+l}) 

N J z 

= e lF{asl) X J — e z +G(logz,a s ,l) ^ (Q ^ 
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where 

oo k 

G(logz,a s ,l) = [logz+l]F{a s [logz+l})-lF{a s l) = 5> fci (logz) V , ( C - 4 ) 

fc=i i=o 

where the coefficients (^(logz) are polynomials in log z. Therefore G represents a 
NLL correction. Expanding its exponential, and integrating term by term in z, we 
get precisely an expression of the form 1 + E(as, I), which is still a NLL correction. 
By taking the derivative of both sides of eq. ( |C.lD with respect to log x we get 



/ diVe -Jviogx + iogiVF( Qs iog7V) =x _E_lQh- x ) e m°st) x [i + E(a s ,l 
2m Jc-ioo ax L 



(C.5) 



which has the form used in section 3. 

In section BL we encounter the Mellin transforms 



/ dNx~ N log m N . (C.6) 
2m J 

We have 

R mk (x) = / x — —, — = — — — T(k,x), (C.7) 

' v ' J 2m N k y ' dk m v ' ; v ; 

where 

™*> = / 5^*" =m Rek>0 - < a8 > 

To prove the above formula, observe that it can always be interpreted according to a 
contour distortion, as the one given in fig. [I]. If k < 1, the contour can be made to 
coincide with the negative real axis, without convergence problems: 



I(k) = -Jdtx 



^ g-fcfOT+logt) _ cc 



2m 



sin(/c7r) 



dte- tlos -r k 

7T 

lX*" 1 sin(A;7r)r(l - k) 
log - 

X J 7T 
fc-1 



1 / 



T(k) 

(using T(l — x)T(x) sin(7rx) = n). This proof is valid for k < 1. On the other hand, 
I(k) can also be given in a form that is manifestly analytical for all Refc > 0, by 
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distorting the contour away from N = 0. By analytic continuation we can therefore 
extend the above result to all the Re k > region. We finally get the result 

= ("J^f e^ 10 ^- 10 ^) . (C.IO) 

From the above formula we also obtain the property 

j 1 dxR m , k (x) = 0, (C.ll) 
J o 

valid for k, m > 0, which follows from 

^og^" 1 (-) dx = F(k). (C.12) 



o \x 
Notice now that 

/ l\ m / l\ m_1 

Rm,i{x) = I -log log- ) -m (-log log- ) 7s + ... . (C.13) 

If we keep only the leading logarithmic term in the expansion, we get 

J R m> i(x) dx & J ^— loglogi^j dx w m! . (C14) 

Observe however that also the integral of the first subleading term is proportional to 
to!, with opposite sign. In fact, if we keep all subleading terms we must get zero, as 
shown earlier. These are precisely the dangerous, factorially growing terms that arise 
when neglecting next-to-leading terms in the x-space resummation formulae. 
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